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Abstract 

We show that the m-dimensional Euler-Manakov top on so*{m) can be 
represented as a Poisson reduction of an integrable Hamiltonian system on a 
symplectic extended Stiefel variety V{k,m), and present its Lax representation 
with a rational parameter. 

We also describe an integrable two-valued symplectic map B on the 4- 
dimensional variety V(2,3). The map admits two different reductions, namely, 
to the Lie group 5*0(3) and to the coalgebra so* (3). 

The first reduction provides a discretization of the motion of the classical 
Euler top in space and has a transparent geometric interpretation, which can 
be regarded as a discrete version of the celebrated Poinsot model of motion and 
which inherits some properties of another discrete system, the elliptic billiard. 

The reduction of B to so* (3) gives a new explicit discretization of the Eu- 
ler top in the angular momentum space, which preserves first integrals of the 
continuous system. 

1 Introduction 

In most publications the integrable m-dimensional Euler top is represented as a flow 
on the cotangent bundle T*SO{m) or on the coalgebra so*{m). 

*AMS Subject Classification 58F07, 70H99, 76B15 
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Recently, an alternative description of this problem as a system on a symplectic 
subvariety of the group product SO{m) x SO{m) was proposed in [1J[3]. 

A first discretization of the free m-dimensional top on T*SO{m) was constructed 
in j24ll2Uj by the method of factorization of matrix polynomials. This discretization 
is represented by a second order Lagrangian correspondence, which does not explic- 
itly involve a time step, it is determined by initial data (a choice of two subsequent 
points on SO{m)). 

On the other hand, in |^ (see also |23]) Suris introduced a concept of an inte- 
grable discretisation of a finite-dimensional Hamiltonian system as a one parameter 
family of integrable Poisson maps parameterized by a time step e, which differ from 
the identity map by 0(e), and whose Poisson structure and the integrals of motion 
differ at most by 0(e) from those of the continuous-time system. 

In the special case where the discretization preserves exactly both the Poisson 
structure and the integrals of motion, one speak of an "exact discretization": one 
has a family of Backlund transformations, which map solutions into solutions and 
are interpolated by a hamiltonian flow generated by some function of the integrals 
of motion of the continuous system. 

A class of implicitly defined Poisson maps so* (3) — > so* (3) discretizing the clas- 
sical Euler top in the space of the angular momentum was indicated in |Bj. The 
maps preserve the energy and momentum integrals of the continuous problem and 
contain explicitly a time step parameter. It was shown that such a map preserves 
the standard Lie-Poisson structure on so* (3) if and only if its restriction onto com- 
plex invariant manifolds, open subsets of elliptic curves, is a shift, which is constant 
on each curve. 

Recently, another integrable discretization of the top on so* (3), which is explicit, 
but does not preserve the integrals of the continuous problem was found in ^1] by 
applying the Hirota method. 

Contents of the paper. Our aim is twofold. First, in Section 2, we propose 
yet another description of the continuous n-dimensional Euler-Manakov top as a 
reduction of a Hamiltonian system on so called extended Stiefel variety V{k,m), a 
symplectic submanifold of dimension km — fc^/2 in M'^"', where 2 < /c < n is an even 
integer. We present a Lax representation of this system with a rational parameter, 
which, in a sense described below, is dual to Manakov's Lax pair found in jl7| . 

The system possesses k/2 commuting symmetry fields TZi generated by Hamilto- 
nians Hi. Its Marsden-Weinstein reduction with respect to the action of the fields 

(k) 

gives rise to a Hamiltonian system on a rank k coadjoint orbit in the coalge- 
bra so*{m), whereas the original Poisson structure in V{k,m) is a pull-back of the 
standard Lie-Poisson structure of so*(m) restricted onto the orbit. The reduced 

(k) 

Hamiltonian system coincides with the Euler-Manakov system on 5^ . In case of 
the maximal rank k, the level variety {Hi = q} C V{k,m) is the group SO{m), 
and the restriction of the original system onto the group yields a flow describing the 
motion of the n-dimensional top in space. 

Second, in Section 3, we present an intertwining relation (discrete Lax pair) 
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generating a explicit A*-depended family of two-valued complex Backlund transfor- 
mations of the variety V(2,3), which preserve the above Poisson structm'e and 
the first integrals of the continuous Hamiltonian system (formula (|3.6j) ). 

The restricton of onto the group 5*0 (3) provides a discretization of the motion 
of the classical Euler top in space and has a transparent geometric interpretation, 
which, in turn, can be regarded as a discrete version of the celebrated Poinsot model 
of motion and which inherits some properties of another discrete integrable system, 
the elliptic billiard (Figure 1). 

On the other hand, the reduction of B'^ onto the coalgebra so* (3) gives a new 
explicit discretization of the classical Euler top, which also preserves its first integrals 
(formula pT^ll '). 

Like the Moser-Veselov correspondence, the both discretizations do not explicitly 
involve a time step and their continuous limits depend on the parameter A*. 

2 Hamiltonian Systems on Extended Stiefel Varieties 
and Rank k Solutions of Frahm— Manakov top 

Recall that the free motion of an m-dimensional rigid body is described by the 
Euler-Frahm equations 

M=[M,n], (2.1) 

where il. £ so{m) is the angular velocity, M G so*{m) the angular momentum of 
the body in the moving frame. Following |18| I21j. these equations are Hamiltonian 
with respect to the degenerate Lie-Poisson bracket on so*{m) 

{Mij,MM}so{m) = SiiMjk - SaMkj + 6kjMji - SikMji (2.2) 

and Qij = dH[M)/dMij. 

The restriction of {•, •}so(m) onto orbits of coadjoint action of SO{m) in so*{m) 
is nondegenerate. A generic orbit parameterized by [m/2] independent Casimir 
functions of the bracket is thus a symplectic variety of dimension m(m — 1)/2 — [m/2]. 

Equations 1)2. 1() are known to be integrable provided M and Q are related as 
[M, a] = [O, b], where a, b are constant commuting matrices, and all the eigenvalues 
of a and b are distinct. The integrability follows from the Lax representation with 
a rational spectral parameter found by Manakov in |17j , or from a hyperelliptic Lax 
pair indicated in [H]. These Lax pairs provide a complete set of integrals of motion, 
whose involutivity can be proved by applying r-matrix theory. 

For the concreteness, in the sequel we consider the case a = diag(ai, . . . ,am), 
b = a?. Then Q = AM + MA, and equations (PTT]) take the form 

M = [M, aM + Ma\ (2.3) 

Apart from this "basic" system, there exists a whole hierarchy of "higher Man- 
akov systems", which are defined by different relations between $7 and M, and which 
commute with (|2.3|) . 
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Below we show that the restrictions of the Frahm-Manakov system on rank k 
orbits of coadjoint representation of SO{m) in so*{m) are closely related to certain 
Hamiltonian dynamical systems on extended Stiefel varieties. Recall that the stan- 
dard Stiefel variety V{k,m) is the variety of ordered sets of k orthogonal vectors in 

(C™) having fixed squares. It is a smooth variety of dimension km — k{k + l)/2 
(see e.g., [7|). 

Namely, as follows from 1)2. the angular momentum in space is a constant 
matrix. Hence, due to the Darboux theorem, in the case rank M = k there exist k 
mutually orthogonal and fixed in space vectors x^^\y^^^ G W^, I = 1, . . . ,k/2 such 
that = ly^'^P = hi and the momentum M can be represented in form 

k/2 

M = Y^ x(') A y(') , that is, M = X'^y - y^ X , (2.4) 
i=\ 

Under the above conditions, the set oikxm matrices Z = (x^^^ y^^^ ■ ■ ■ x^^^"^^ y{'^/2))T 
forms the extended Stiefel variety V{k,m). In contrast to the standard Stiefel 
variety, for each index Z, the absolute values \x'^^'^\ = are not fixed. Thus, 

V{k, m) is of dimension km — A:^/2, and the kxm components of Z play the role of 
excessive coordinates on it. 
Let 

k m 

w = tr {dX A dy'^) = YY^ dxf A dyf^ 

1=1 i=l 

be the canonical symplectic structure on the space M'^'" = x^^\y^^\ • • • , x^^/'^\y'^^f'^^) 
and let a) denote the restriction of 2-form uj onto V(/c, m) C R'''". The latter subva- 
riety is defined by conditions 

- = 0, (x«,x(^)) = = 0, yX^ = Q, (2.5) 

/,s = l,...,A;/2, 

which consist of k"^ /2 independent scalar equations fs{x,y) = 0. The matrix of 
standard Poisson brackets of the constraint functions fs in M'^™ is nondegenerate. 
It follows that 2-form O is also nondegenerate and the extended Stiefel variety is 
symplectic. 

Since the vectors are fixed in space, in the frame attached to the top they satisfy 
the Poisson equations x^^-* = —Qx^^\ y^^^ = —fly^^\ Q G so{m), which imply 

x = xn, y = yn. (2.6) 

As above, we put O = aM + Ma, a = diag(ai, . . . ,Om) and define a dynamical 
system on V{k,m), which is generated by (|2.6j) . (|2.4|) : 

X = x[a{x^y -y^x) + x^ya], 
y = y[aix^y - y^x) - y^Xa] . 
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Theorem 2.1 1) Under the substitution ^2.4\l solutions of the system \2. 1\ give 
rank k solutions of the Frahm-Manakov system i2.cl\) on so* (m) . 

2) Up to the action of the discrete group generated by reflections (X,y) {—^, —y), 
the system \2. ?D is described by k x k Lax pair with rational parameter A 

L{X) = [L{X),A{X)], L,Aespik/2), AgC, (2.8) 

i=l 



-lyT 


-X{X1- 




-lyT 


y{xi- 


















\-yiyJ 





(X{a + XI)y'" X{a + XI)X^\ 

where Xi = (^x\^\ . . . , xf'^'^^^ (respectively yi = . . . , x-*^''^^^ ) is i-th 

column of X {respectively ofy), and I is the unit m x m matrix. 

Proof. The first statement follows directly from the derivation of the system (|2.7|l . 
Further, we calculate the derivative L{X) by virtue of equations ()2.7() . In view of 
matrix relations (AI — a)~^a = A(AI — a)~^ — I and yX^ = Xy^ = 0, the derivative 
coincides with the commutator in ()2.8|) . q 

Remark 2.1. Notice that the entries of matrices 

$(A)L(A), $(A)A(A), where $(A) = (A - ai) • • • (A - a„) 

are polynomials in A, and, under the substitution (|2.4I) . the coefficients of the char- 
acteristic polynomial \^{X)L{X) — wl\ can be expressed in terms of Mij only as 
follows 

k _ 
|'«;I-L(A)| =u;'= + ^u;^-^^''-^(A)J/(A,M), I = 2,4:, . . . ,k, 

1=2 

%{X, M)=Y, TT T^^T 7 \M\\ , (2.11) 

^ {X - ai^) ■ ■ ■ {X - a^^) 

where \M\\ are diagonal minors of order / corresponding to multi-indices I = 
{ii,...ik} C {l,...,m}, 1 < ii < • • • < < m. Notice that the leading coeffi- 
cients Hi fn-iiM) = form a complete set of Casimir functions on so*{m). 

The k X k matrix L(X) in (|2.9j) belongs to a wide class of Lax operators of the 
form 

X-ai 

1=1 
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where Y G gl{k) is a constant matrix and Gi,Fi are ki x k matrices. Such Lax 
matrices can be regarded as images of moment maps to the loop algebra gl{k), and 
integrable systems generated by them have been studied in the series of papers [2 
mElE] in connection with the duality to so called rank k perturbations of constant 
diagonal matrices of dimension nxn (following Moser [Hi)- In particular, the kx k 
Lax matrix H2.9|) is dual to the nxn Lax matrix in the Manakov representation, 

£{^i) = a + -{X^y-y^X) = a + -M, (2.12) 

in the sense that under the relation ()2.4|) the spectral curves \L{X) — = and 
|>C(//) — wl\ = are birationally equivalent and the parameter A plays the role of 
the eigenvalue parameter for (|2.12l) . The characteristic polynomials of the dual Lax 
matrices are related by the Weinstein-Aronjan formula (see [Tj). 

Remark 2.2. The matrix ^(A) in (|2.10l) can be represented in form 

A(A) =[A— +2cj,(A)L(A)]+ + (ai + --- + aOLo, Lo = ^_yyT q 

where [ ]+ denotes the polynomial part in A of the expression. Notice that the Lax 
equation L = [L, Lq] describes the vector flow 

= y(0 = -(y(0^y(0)a;{0^ l = l,...,k/2. (2.13) 

For each index TZi, equations (|2.13|1 generate rotations TZi in 2-planes spanned by 
the vectors x^^\y^^\ which leave the momentum M invariant. 



Let {•, •} be the Poisson bracket on V{k,m) obtained as the Dirac restriction of 
the standard bracket in M'^™. Symplectic properties of our system are descibed by 

Proposition 2.2 The dynamical system \2. 7| ) on V(A;, m) is Hamiltonian with re- 
spect to {•,•} with the Hamilton function H{X,y) = — |tr(M^(Af, 3^)^). In the 
abundant coordinates X,y it admits the canonical representation 



dH 



Ji) 

~ - ,(0 



dyl 



.{I) 



dH 



V{k,m) " dxf^ 



(2.14) 

V{k,m) 



i = 1, . . . , n, 1 = 1,..., k/2. 

Proof. The equivalence of equations 1)2. 14() and ()2.7|) on V{k,m) is verified by direct 
calculations. Next, according to the Dirac formalizm, the standard bracket and {•, •} 
are different by terms containing {fs, H}. The latter equal zero since the constraint 
functions fs given by ()2.5|) are invariants of the flow generated by H{X, y) on M'^'". 
Hence, equations (|2.7|) or (|2.14|) are Hamiltonian with respect to {•, •}. q 

Rotations IZi given by (|2.13j) are generated by the Hamiltonians Hi{x,y), the 
restrictions of the functions ^{x^^\x^''^){y^^\y^^^) onV{k,m). Clearly, these functions 
are first integrals of the system 1)2. 7|1 and moreover they commute with H. 
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Let us fix the values of the Hamiltonians by putting 



(x(^), x(') ) = = /i/, /i/ = const /O, / = l,...,A;/2. 

These conditions define the customary Stiefel variety V{k,m). Under the substitu- 
tion ()2.4() . the factor variety V{k,m) /{TZi, . . . ■,'R-k/2} coincides with a rank k coad- 

joint orbit C so*{m) of dimension k [m — ^ — k, which is parameterized by 
the constants hi, ... , Notice that = + • ■ ■ + 

(k) 



Theorem 2.3 



1) Under the map V(A;,m) 



so*{m) is the push-forward of the bracket {•,•}. 



Sf^' , the Lie-Poisson bracket on 



2) The Poisson (Marsden-Weinstein) reduction of the system \2. 7| j obtained by 
fixing values of Hi{x,y) and by factorization bylZi, I = l,...,fc/2 coincides 
with the restriction of the Frahm-Manakov system with Hamiltonian 
H{M) = \ J2i<j{ai + aj)M^j onto the orbit 5^ \ 

Proof 1). In view of ^l^i, (1^ . 

{Mi,{x,y),Mki{x,y)} 



{M,j,MM}so{n){^,y), 



i.e., the canonical bracket on M'^™ is the pull-back of the bracket {•, •} 
so*{m). On the other hand, on V{k,m), 



so{n) on 6;^ C 



{Mi,{x,y),Mki{x,y)} = {Mij{x,y),Mki{x,y)}, 

since for any i,j,s, {Mij{X ,y), fs{X ,y)} = 0. This proves item 1). 

2). By item 1) and Proposition 12. 2( the Poisson reduction of system (|2.7|) onto 
is described by the Lie-Poisson bracket {•, ■}so{n) ^^'^ the Hamiltonian H{M) = 
H{X,y) = X^j<j(aj + aj)Mfj, i.e., it is the corresponding restriction of the Frahm- 
Manakov system, q 

(k) 

The reduced system on the orbit 5^ is integrable and its generic invariant 
manifolds are tori of dimension idim5^'^'* (see, e.g., HE])- On the other hand, 

the preimage of a generic point M £ 5f ^ in V{k ,m) is a /c/2-fold product of circles 
5^ X • • • X 5^ (in the complex case C* x • • • x C*). This implies that the original system 
on V{k, m) has generic invariant tori of dimension ^dimcS^*^^ + k/2 = (m — k/2)k/2, 
i.e., a half of dimension of the symplectic manifold V{k,m). Hence, the original 
system 1)2. 7|) is also integrable. 

To get a global view on the above manifolds, we represent them in the following 
diagram, with the dimension indicated above, where arrows denote the correspond- 
ing relations (restrictions or factorizations). 



fs=0 



V{k, m) 



;(0|2 = |j,{i)|2=;,j 



V{k, m) 



n 



S 



(k) 



k X m 











k[rn--j 





k { m 

2 



k [ m 

2 
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Remark 2.3. In the case of maximal rank k {k = m or k = m — 1), when is 
a generic coadjoint orbit Sh, the Stiefel variety V(/c, m) is isomorphic to the group 
SO{m). Then the following commutative diagram holds 



V(k, m) 



7^ 



so*{m) 



|xW|2 = |y(0|2=/i; 



Hl,m-liM) = Cl 



SO{m) 



n 



where the values {q} of nonzero Casimir functions Hi^^_i{M) correspond to the con- 
stants {/i;}. The mapping 5(5 ^ > 5^ can be regarded as a multi-dimensional 

analog of the Hopf fibration S0(^) > S"^- '^^^ restriction of the system (j2.7|) 

onto V(fc,m) yields an integrable flow on the group SO{m) which describes the 
motion of the Frahm-Manakov top in space for the chosen angular momentum. 

For m = 3 such a flow was considered in jl5l I16j from the point of view of its 
hydrodynamical interpretation. 



A generalization of the Chasles theorem. If the rank k is not maximal, then 
the components oi X ,y themselves are not sufficient to form a complete set of 
coordinates on SO{m) and to determine the position of the top in space uniquely. 
However, in this case one can make use of the following geometric property described 
in [H] . Let us fix a part of constants of motion by putting in l|2.11j) 

lk{s,M) = co{s - ci)---{s - Cm-k), Co, ci, ... ,Cm_fc = const (2.15) 
and consider family of confocal cones in R™ = {Xi, . . . , X„) 

Q[s) = \^^ + ... + ^^ = q\. (2.16) 
[ s - ai s- an J 

Let A C be a /c-plane spanned by the orthogonal vectors x^^^ , y^^^ , • • • , x^^^"^^ , y^^/"^^ . 

Proposition 2.4 (0). 

1) . Under the motion of the Frahm-Manakov top with constants \2.15\) the k-plane 

A is tangent to the fixed cones Q(ci), . . . ,Q{cm-k)- 

2) . Let (p^^^ be a normal vector of the cone Q{ca) at a point of the contact line An 

Q{ca). Then the vectors . . . , (f)(^-^) together with y^-^^, . . . , x^^/'^\y''^/'^^ 
form an orthogonal frame in M™ which is fixed in space. 

For fixed polynomial Zi^{s,M), the vectors (f)^''^ can be calculated in terms of 
x^^^ , y^^^ and, thereby, the position of the top in space is completely determined. 
Proposition 12.41 defines a single- valued map V{k,m) — > SO{m) under which generic 
invariant tori of dimension (m — k/2)k/2 on V{k,m) become tori of the same di- 
mension on the group SO{m). 

Note that the above proposition generalizes the celebrated Chasles theorem on 
the propery of the tangent line to a geodesic on a quadric. 
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The rank 2 case. In the simplest case k = 2 the angular momentum can be 
represented in form 

M = xAy, X = X^^^ = {xi,...,Xmf , y = yW = {y^, . . . ,y^f (2.17) 

and equations H2.7() describe a Hamiltonian system on the extended Stiefel variety 
V(2,m) = \{x,y) 



\x\ = \y\, {x,y) = 

X = -{y,ax)x + {x,ax)y + ay{x,x), 
y = -{y, ay)x + (x, ay)y - ax{y, y) 

with the Hamiltonian 

1 1 1 

H = -{x,ax){y,y) - {ax,y){x,y) + -{y,ay){x,x) = - ^(oi + aj)M^-. 



Equivalently, this system describes the evolution of fixed orthogonal vectors x, y 
in a frame attached to the m-dimensional body. The system admits the following 
2x2 Lax pair arising from (|2.8|) . 



L(A) = [L(A),^(A)], (2.19) 



L(A) = $(A) 
-4(A) = 



i=l X—ai Z-ti=l \—ai 

Em Vj sr^m Xj-yj 

i=l X—ai Z^i=l X—ai 



YJILii^ + 0'i)xiyi -YlT=i{>^ + ai)x\ 
YT=i + (^i)y'i YT=i (-^ + o-i)xiyi 

^>(A) = (A-ai)---(A-a„), 

The Lax representation (|2.19j) was first indicated in [2j, where it was shown to 
be dual to an n x n Lax pair for the rank 2 case found by Moser in |19j . 

In view of relation (|2.17() . the characteristic polynomial \L{X) — ^I| for 1)2. 19(1 
can be written in form $(A)T2(A, M) + /i^, where X2(A, M) is the family of quadratic 
integrals defined in 1)2. 11() . 



T2(A, M) = V ^[^] -Ml 

f^j (A-a,)(A-a,) 

= A'"-2//2,m-2(M) + X^~^H2,rn-3{M) + • • • + H2o{M). (2.20) 



r2 



Notice that H2^m-2iM) = Yli<j ^ij — {y^y){x,x) is a Casimir function of the 
standard Lie-Poisson bracket on so* (m) . With respect to the Poisson bracket { , } 
on V(2,m), this function generates permanent rotations of the top in the fixed 2- 
plane A = span(y,x), which leave the components of M invariant. 
Let us fix the constants of motion by putting 

X2(A,M) = co(A - ci) • • • (A - Cm_2), Co, ci, ... ,Cm-2 = const. (2.21) 
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This defines hyperelliptic spectral curve in = (A, fi) of genus g = m — 2 



C = {/.2^-co$(A) (A-ci)---(A 



2)}- 



(2.22) 



As noticed in the real generic ((7 + l)-dimensional invariant tori of the system 
can be extended to open subsets of generalized Jacobian varieties Jac(C, oo±), which 
are extensions of the customary (7-dimensional Jacobian Jac(C) by C* and which can 
be regarded as the factor of by the lattice generated by {2g + 1) independent pe- 
riod vectors of g holomorphic differentials Oi . . . ,LDg and a meromorphic differential 
of the third kind f^oo± having a pair of simple poles at the infinite points oo± on 
the curve C. 

The coefficients of the matrix polynomial L{X) are meromorphic functions on 
Jac(C, 00-1-), whereas the components of the momentum Mij and the normal vectors 
are meromorphic on a covering of the Jacobian Jac(C) itself (the C*-extension 
is factored out by the action of 7^ = 50(2)). 

In the classical case m = 3 the curves C become elliptic ones and generic in- 
variant tori in V(2, 3) and in SO{3) are 2-dimensional. An explicit solution for the 
components of the rotation matrix in terms of theta-functions and exponents was 
first given in JSI (see also ^2^). Since now rank M =2 in the generic case, the above 
commutative diagram takes the form 



V(2,3) SO{3) 



SO{2) 



SO* (3) 



SO{2) 



{M,M)=h'^ 



02 



S\ being the coadjoint orbit (2-dimensional sphere) corresponding to the constant 
h. 



3 Backlund transformation on V(2,3), 50(3) 
and discretization of the clasical Euler top 

A first integrable discretization of the m-dimensional Euler-Manakov top was con- 
structed in |24| I2()j by the method of factorization of matrix polynomials. It was 
represented by the correspondence {yt,,M) (Q,M), ft £ SO{m),M E so*{m), 
which, in our notation reads 

M = n'^A- An, M = nMn'^. (3.1) 

Given M, the new matrix Q is found from equation M = Cl^'^A — ACi, whose solution 
is not unique. 

In given section we describe a symplectic map 8x^ : V(2, 3) V(2, 3), 8\^{x, y) = 
(x, y) governed by an arbitrary parameter A* G C, which preserves the first integrals 
of the continuous system 1)2. 18() and whose restriction to each generic complex torus, 
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generalized Jacobian Jac(C, oo±), is given by shift by the 2-dimensional vector 



S= — , , i?± = (A*,±/i*), 

E± being involutive points on the elhptic spectral curve, the simplest case of 1)2. 22p . 

C = {i? = -Co (A - ai)(A - a2)(A - 03) (A - ci)}. 

Note that 5" is a correctly defined vector in the generalized Jacobian: under a change 
of integration path on C it increases by a period vector of Jac(C,oo±). We also 
emphasize that here A* is a constant parameter, whereas the conjugated coordinate 
/i* depends on the equation of the curve. 

Backlund transformation on V(2, 3). As shown in ^Uj by applying an addition 
theorem for a class of meromorphic functions on generalized hyperelliptic Jacobians, 
such a map admits intertwining relation (discrete Lax pair) 

Z(A)M(A|A*) = Af(A|A*)L(A), (3.2) 

M(.,A-,^(-"<\-;)-V„„_V,-,)> 

B(\*\ = + _ + a\alal{x,{a-X*l)-^y) 

^ ' LuiX*) ala*^a*^{x,{a-X*I)-^x) ' 

1 O-i — O-i ~ ^ J 



a 



dX 



x=x* 



where L{X) is defined in ()2.19|) and L(A) depends on the new variables x, y in the 
same way as L{X) depends on x, y. In view of (|2.2()j) for m = 3, 



HX*)^al{xiyj-Xjyi)\ (i, j. A:) = (1, 2, 3). (3.4) 



k=l 



Now putting in 1)3. 2() subsequently A = 01,03,03 and calculating the matrices 
M(oj)-L(oj|x, y)M^^(oj), we find 



xl 



{yi + I3xi + a{ai - X*)xi)'^ 



a2(oi - A*)2 
S^Xj — a(oj 

a2(oj — x*y 



Vi = T^T- TIVo ' (3-5) 



~ ~ _ {yi + I^Xj + a{ai - X*)xi) iPyi + P'^Xj - a(oi - X*)yi) 

a^(oj — A*j^ 

From here the new variables can be recovered up to the action of the group generated 
by refiections {xi,yi) — > {—Xi,—yi). Imposing the condition of the existence of a 
continuous limit (see below), we choose the following relations 

~,_ JJi + (^ ^r._y.-_^k±M i-12 3 (3 6) 
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These expressions together with p.3p . (|3.4|) describe the map Bx^ : V(2, 3) V(2, 3) 
in an explicit form. Since a generic parameter A* corresponds to two values of 
the map is generally two-valued. 

Geometric model. The restriction of the map onto the group 5*0 (3) admits a 
transparent geometric interpretation, which can be regarded as a "discrete version" 
of the kinematic Poinsot model (see, e.g., ^2S]). Namely, let \x\ = \y\ = 1 and let 

R=\\xy xAy\\£ SO{3) 

be rotation matrix defining a position of a rigid body in space. We attach to the 
body a cone K2 = {{X, (a — A*I)~^X) = 0}, which is fixed in the body frame 
{Xi, X2, Xs), and assume that 

< ai < 02 < 03, ai < A* < a2 or 02 < A* < 03. (3.7) 

Under these conditions the cone is real and regular. Let 11 be 2-plane spanned by 
X, y, which is thus fixed in space and orthogonal to the momentum vector M = xf\y. 
Assume also that x, y are such that 11 has a nonempty real intersection with the 
cone K2 along lines Li, L2. One can show that under this condition the coordinates 
/i* defined in (|3.4j) and the parameters a, /? are real. 

Theorem 3.1 Let Ki = {{X,{a - hl)-^X) = 0}, /i = const he a unique cone 
attached to the body such that it is confocal to K2 and tangent to the fixed plane 11. 
Then then new position of the body defined by the rotation matrix R = \ \xy x f\y\\ 
and expressions \S. 6)) is obtained from the original position by rotating the cones 
Ki,K2 about axis Li or L2 until Ki again touches 11. 

This geometric construction is illustrated on Figure 1. In the new position R 
determined by rotation about L2, the cone K2 intersects 11 along L2 and another 
line -L3. Then the next iteration is generated by rotation about L2 or L3. The 
two-valuedness of the map Bx^ is now related to the possibility of rotation about 
two different axes in M^. It follows that iV-th iteration of the map is only (A^ + 1)- 
valued, not 2^-valued. By fixing a sign of /i* in (|3.4|) . Bx* becomes single- valued 
and generates a sequence of points on SO {3). 

The geometric model was first proposed in 10, as a certain limit of a kinematical 
model of motion of 4-dimensional Frahm-Manakov top in space. 
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Figure 1 

Remark 3.1. As follows from the map ^B^* admits particular solutions, for 

which the vector M = x A y remains to be an eigenvector of the inertia tensor 
A = diag(ai, 02, 03), whereas x,y themselves rotate by a fixed angle in the plane A. 
Such solutions can be regarded as analogs of stationary permanent rotations of the 
classical Euler top about its principal intertia axes. 

Continuous limit. Note that when x, y are chosen such that 11 is (almost) tangent 
to the cone K2 (c A*), Ki and K2 confluent and, according to the above model, 
the cone Ki = K2 is rolling without sliding over the fixed plane thus giving a 
continuous limit motion on the group 50(3). 

From the algebraic geometrical point of view, in the above limit the points 
on the spectral curve C come together to a branch point Eq = (c, 0) and the 
shift vector S on the generalized Jacobian tends to zero. Let e be a small complex 
parameter. Setting /u* = e, A — A* = const • in (|3.3|1 . we have the expansions 

{x,a-'y) , ^, , _ 1 , _i « , _ dfi\X) 



/3 = - ' +0(e), a = -x{x,a-^x) + 0{1 



\=\* 



{x,a ^x) ' e ' ' d\ 

where x is a real nonzero constant. Now we set 

X = X + xe + O(e^), y = y + ije + 0{e^). 

Substituting the above expansions into (|3.6j) . then matching the coefficients at e in 
both sides and taking into account relation 

(x, (a*)-iy)2 = (x, (a*)-ix)(y, {a*)-\) + 0{e^), 

we obtain the following differential equations describing the limit flow on a subset 
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of V(2,3) 



X 



y 



deta* 
deta* 



{x, {a*r'y){a*)-'x - {x, {a*)-'x){a*)-'y 



1 



1 



X Aa*{x A y), 



(3.8) 



= -yAa*{xAy), 



a* = a- X*I. 

These equations are Hamiltonian with the Hamilton function 
det a 



H 



2k 



{x,{a*r'x){y,{a*r'y) - {x,{a*)-'yf 



i=l 



Notice that this function equals zero on the limit continuous flow. The restriction 
of this flow on Jac(C, oo±) is tangent to the curve C C Jac(C, oo±) at the point Eq. 

The above asymptotic of a, (5 explains the specific choice of sign of x, y made in 
the passage from relations (|3.5() to the map ()3.6() . 

Proof of Theorem Vd.ll The condition for 11 to be tangent to the cone 
Ki = {{X, (a - hI)-'^X) = 0} has the form 

3 

^(/i - ak){xiyj - xjyif = l2{h,M) = 0. 
fc=i 

Comparing this with the family of integral 1)2. 21() for m = 3, we conclude that 
h = ci, which is constant under the map, hence the plane spanned by x,y is again 
tangent to Ki. 

Next, any translation in 5*0 (3) is represented as a finite rotation about an axis 
in M^. As follows from relations 1)3. 6|) . y + (3x = y + fix, hence the line along the 
vector i = y + /3x is invariant of the action of Bx^ on and therefore represents 
the axis of such a rotation. Finally, in view of H3.3|) . we have 



{(i,{a*)-H) = {y,{a*)-'y) 



+ 



Li2{X* 

(/x*)2 + 2/i*Ln(A*) + Lfi(A 



{x,{a*)-'x) 



1 



Lu{X 



which equals zero by virtue of H3.4() . Hence {i, (a*) = 0, which imply that the 
vector i lies on the cone i^2- This establishes the theorem, q 

Remark 3.2. When the attached cone K2 does not have real intersection with 
H = span(a;,y), the coordinate /i* is imaginary and, according to 1)3. 6() . (|3.3|) . the 
new values x, y are complex. As a result, under the reality conditions (|3.7|) the map 
is real only on the subset ^ C V(2, 3) defined by unequality 

3 

^(A* - ak){xiyj - Xjyif < 

k=l 

On the boundary of fH, the map tends to the identical one. 



14 



Reduction to the coalgebra so* (3). Under the factorization by rotations of 
TZ = SO{2), the transformation 8x^ induces a map Bx^ : so{3)* so{3)* such that 

M = Bx*M{x,y) =xAy. 

The latter map is correctly defined, i.e., it does not depend on a concrete choice 
of vectors x, y giving the same M. It preserves the first integrals of the classical 
Euler top on so* (3) and its generic invariant manifolds are open subsets of 4- fold 
unramified coverings of the complex torus Jac(C) = C. The restriction of Bx* onto 
Jac(C) is given by shift by the holomorphic integral e = J^'^ dX/n, which thus 
depends only on the constants cq,ci. According to a theorem in [S], this implies 
that the map iS^* preserves the standard Lie-Poisson structure on so* (3). 

Proposition 3.2 Vectors M, M satisfy the following symmetric relations 

M - M = x{M + M) Aa(M + M), (3.9) 

/ — — Jl- {M,M)/co 

X = Y 2(oM + aM, aM + aM) ^ ^ == , (3.10) 

^l + (M,M)/co 
(M, a*M) = (M, a*M) = -{M, a*M), (3.11) 

where, as above, cq = {M,M) = {M,M). 

Relation 1)3. 9|) was previously obtained by another method in [S], as an implicit 
map describing a Poisson discretization of the Euler top in so* (3). 

Proof of Proposition 13.21 In view of relations (|3.6|) , we find 



X Ay = M 



xAy + -(a*y^£Ae, (3.12) 
a 



where, as above, I = y + f3x, a* = a — A*I. Note that vector (a*)~^ i is normal to the 
cone K2 at a point of the intersection line L2 or Li. Hence, M — M is orthogonal 
to {a*)~^ i and i. Next, since i lies in the planes 11,11, this vector is orthogonal to 
M,M. This, together with the equality \M\ = \M\ implies that the sum M + Af is 
parallel to (a*)"^ i and a*(M + M) is parallel to i. As a result, (|TT^ implies 

To find factor x, we first introduce angle (j) between vectors M and M. Since 
\M\ = \M\, the vectors M — M and M + M are orthogonal, and we have 



~ 1~ 6 1— \^l-{M,M)/co 

\M -M\ = -|M + M|tan- = -\M + M\ ^ ^ — 

2 2 ^ ^l + (M,M)/co 

On the other hand, since a*(M + M) is orthogonal to M,M, from 1)3. 9|1 and the 
properties of the vector product we deduce 

\M - M\ = x\M + M\ \a*{M + M)\. 
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Comparing the right hand sides of the above two relations, we obtain ()3.1U() . 

The first equahty in (|3.11|) holds because the map Bx^ preserves the first integrals 
of the Euler top. Next, since the vector a*{M + M) lies on the cone K2, we have 
{{M + M),a*{M + M)) =0. Expanding this and using the first equality in (fXTTj) 
yields the second equality. □ 

Remark 3.3. The fact that the difference M — M is orthogonal to (a*)^^ £ and 
i implies that the angle between the normal vector (a*)^^ £ and the plane H equals 
the angle between (a*)~^£ and H. This property can be regarded as a projective 
version of the Birkhoff condition of elastic impacts, hence the geometric construction 
of Theorem 13.11 illustrated in Figure 1 describes a projective analog of the plane 
elliptic billiard. (Note that no any plane section of the cones and of the sequence of 
n gives such a plane billiard.) 

To obtain the map ^a* : so* (3) — > so* (3) in an explicit form, we use the fact 
that the vector i satisfies the system of homogeneous equations 

{i,{a*)-^£) = 0, {£,M)=0. 

One of its solutions, i = (^1,^25^3)"^) normalized by the condition £3 = 1, has the 
form 

- , MiMs-^/Da*M2 „- , M2M3 + /PatMi - 

ti = — ai T T — , to = —cin o o — , = 1 (3.13) 

^ alMf + a*2Ml ' ^ a^M^ + a*2Ml ' ^ ^ ^ 

D = fr^; 3.14 

Substituting this instead of £ into (|3.12j) and symmeterising obtained expressions, 
we arrive at the following relations 



Ml 


- Ml 


= X(«2 


M2 


-M2 


= X(«3 


Ms 


-Ms 





as) (Ai M2M3 - {a\fMl{alMl - alMl) \ 
ai) (a2 Ml Ms - Vd {alfMi{alMl - ajMf ) ) , (3.15) 



where 

A. = n*n* 



a*alM^Ml - {alfuf, k) = (1, 2, 3) 



and X is £^ common factor. 

Next, multiplying the both sides of ()3.15() by OiMi, a2M2, OsMs respectively, 
then summing and using the second equality in (|3.11|) . we find the factor x iii form 

X = r T ' (3-16) 



Y,ai{a2 - as) 



A1M1M2M3 - ^/D{a\fMf{alMl - a^Ml) 



where the summation in the denominator ranges over the three terms obtained by 
the cyclic permutations of indices (1,2,3). Thus the right hand sides of (|3.15j) are 
homogeneous expressions of degree -1 in Mj. 

We summarize our results on the map B\^: in the following theorem. 
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Theorem 3.3 The map Bx* : so* (3) so* (3) given by /TO^I) . ^TTB, 

preserves the first integrals of the continuous Euler top, as well as the standard Lie- 



Its restriction onto each regular elliptic curve C is represented by the shift by the 
holomorphic integral e = jj'^ dX/fi. Near the boundary ofD\ the energy constant ci 
tends to A* and the shift e tends to zero. The continuous limit of B\^, coincides with 
the Euler equations M = [M, oM] . 

Remark 3.4. Since generic solutions Mi(t), M2(t), M3(t) of the classical Euler 
top are proportional to the elliptic functions sn(u), cn(u), dn(ii), with u = const • t 
and the map B\^, results in adding e to the argument u, the relations H3.15() . H3.14() . 
(|3.16|) can be regarded as a set of explicit addition formulae for these functions. 
However, in this case neither the moduli of the curve C nor the shift parameter are 
fixed: they depend on initial position of M in so* (3). 
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